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ESTIMATING  THE  ROTATION  RATE  OF  MESOSCALE  OCEAN  RINGS 

FROM  SATELLITE  ALTIMETRY 


INTRODUCTION 

The  U.S.  Navy’s  thermal  analysis  system  assimilates  data  into  a  three-dimensional  depiction  of 
the  ocean  thermal  structure.  The  input  data  for  the  analysis  includes  bathythermographs,  ship  and 
buoy  temperatures,  and  satellite-derived  sea  surface  temperatures.  In  addition,  a  key  input  to  the 
system  is  a  map  of  the  surface  locations  of  fronts  and  rings,  together  with  estimates  of  such 
parameters  as  the  ring  rotation  rate.  This  map  is  assimilated  into  the  analysis  via  feature  models. 
A  feature  model  is  a  schematic,  diagnostic  model  that  depicts  the  typical  three-dimensional  thermal 
structure  of  fronts  and  rings,  and  is  used  to  infer  the  subsurface  thermal  structure  from  the  surface 
location  of  fronts  and  rings.  It  is  a  powerful  tool  for  supplementing  sparse  data  with  a  knowledge 
of  the  oceanography  of  mesoscale  ocean  features. 

Ring  parameters  are  used  to  tune  the  ring  feature  model  to  agree  with  the  characteristics  of  a 
particular  ring.  Experience  has  shown  that  the  ring  rotation  rate  is  a  key  parameter  for  tuning  the 
ring  feature  model.  This  report  presents  a  technique  for  using  satellite  altimetry  data  to  estimate 
the  ring  rotation  rate.  The  technique  is  applicable  even  if  the  altimeter  track  does  not  cross  the  ring 
center.  Applicability  to  tracks  that  do  not  cross  the  ring  center  is  important  because  the  study  of 
Tournadre  (1990)  suggests  a  less  than  5%  probability  that  a  given  ring’s  center  will  lie  on  a 
GEOSAT  or  TOPEX/POSEIDON  altimeter  ground  track.  He  also  considers  the  case  of  an  altimeter 
crossing  within  half  a  radius  of  the  ring  center.  The  percentage  of  Gulf  Stream  rings  that  have  a 
probability  greater  than  50%  of  being  detected  (distance  between  the  ring  center  and  altimeter  track 
is  less  than  half  a  ring  radius)  is  greater  than  83%  for  the  warm-core  rings  and  45%  for  the  cold- 
core  rings. 

Section  2  presents  some  background  on  thermal  analyses,  front  and  ring  maps,  and  feature 
models.  The  technique  for  estimating  the  ring  rotation  rate  is  discussed  in  Sec.  3.  Examples  of  the 
application  of  this  approach  are  presented  in  Sec.  4.  Section  5  discusses  how  well  the  estimates 
agree  with  rates  that  have  been  reported  in  the  literature,  as  well  as  the  potential  for  assimilating 
drifting  buoy  data  into  a  composite  estimate.  Conclusions  are  set  forth  in  Sec.  6,  and  Sec.  7  has 
recommendations  for  future  work. 


BACKGROUND 

Two  key  aspects  of  monitoring  the  ocean  thermal  structure  are  data  and  a  thermal  analysis 
system.  The  data  are  from  disparate  sources  irregularly  sampled  in  space  and  time,  each  with  its 
own  characteristic  quality.  The  thermal  analysis  system  assimilates  these  data  into  a  regularly 
gridded  map  of  the  ocean  thermal  structure. 
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The  U.S.  Navy  uses  optimum-interpolation-based  (Gandin  1965;  Bretherton  et  al.  1976;  While 
and  Bernstein  1979)  thermal  analysis  systems  for  making  three-dimensional  analyses  of  the  ocean 
thermal  structure.  One  system  is  known  as  the  Optimum  Thermal  Interpolation  System  (OTIS) 
(Clancy  1987;  Clancy  et  al.  1990).  This  sy.stem  assimilates  bathythermograph  (BT)  data,  front  and 
eddy  maps,  and  sea  surface  temperature  (SST)  data  from  satellites,  ships,  and  buoys.  In  addition, 
an  OTIS  analysis  is  used  to  initialize  a  forecast  of  the  upper  (400  m)  ocean  thermal  structure  by 
the  Thermodynamic  Ocean  Prediction  System  (TOPS)  (Clancy  and  Poliak  1983).  The  thermal  analysis 
can  also  be  used  to  initialize  a  dynamical  ocean  forecast  model  (Robinson  and  Walstad  1987;  Fox 
et  al.  1993).  OTIS  is  presently  operational  in  some  regions  at  the  Navy's  Fleet  Numerical  Ocean¬ 
ography  Center  (FlsJOC).  The  Modular  Ocean  Data  Assimilation  Sy.stem  (MODAS)  is  a  more  modular 
and  user-friendly  version  of  OTIS  that  was  developed  at  the  Naval  Research  Laboratory.  MODAS 
is  being  implemented  in  several  areas  for  use  at  the  Naval  Oceanographic  Office  (NAVOCEANO). 

A  major  impediment  to  the  effective  use  of  any  thermal  analysis  system  is  the  sparsity  of  data 
available  for  monitoring  the  mesoscale  ocean,  especially  at  depth.  The  BT  is  the  principal  source 
of  subsurface  temperature  data,  but  there  are  only  about  200  BTs  per  day  for  the  world's  oceans 
(Clancy  1987).  SST  data  are  obtained  from  ships  and  buoys.  Although  ship  and  buoy  SST  data  are 
more  abundant  than  BT  data,  ship  data  have  larger  errors.  Satellite  data  offer  the  prospect  of  a  cost- 
effective  system  for  the  long-term  monitoring  of  the  ocean  environment  (Hawkins  et  al.  1989). 
Satellite-derived  SSTs  have  been  used  routinely  for  a  number  of  years  by  the  U.S.  Navy  as  data  for 
OTIS.  These  data  are  unobtainable,  however,  in  cloud-covered  areas. 

More  recently,  mesocale  fronts  and  rings  have  been  located  using  satellite  altimeter  data  (Lybanon 
et  al.  1990;  Cheney  and  Marsh  1981).  Altimeter  data  are  particularly  useful  for  locating  submerged 
cold-core  rings  that  have  no  SST  signature  (Lybanon  et  al.  1990).  However,  the  width  of  the 
altimetric  footprint  is  very  narrow.  When  the  satellite  is  in  an  exact  repeat  mode,  the  ensemble  of 
groundtracks  form  a  grid  that  is  periodically  repeated.  The  spacing  between  the  tracks  can  be 
reduced  only  at  the  expense  of  increasing  the  time  between  repetitions  of  the  tracks.  In  the  case  of 
the  GEOSAT  Exact  Repeat  Mission,  which  was  optimized  for  oceanographic  applications,  the  grid 
has  a  spacing  of  1.4754  degrees  at  the  equator  and  a  17-day  repeat  cycle  (Born  et  al.  1987).  Hence, 
several  days  can  elapse  before  a  ring  is  revisited  by  the  altimeter. 

An  important  input  to  OTIS  is  a  map,  known  as  the  bogus,  of  the  surface  location  of  mesoscale 
ocean  fronts  and  rings.  This  map,  together  with  estimates  of  various  ring  parameters  and  front 
characteristics,  is  produced  by  NAVOCEANO.  The  bogus  is  produced  primarily  by  examining  the 
thermal  gradients  found  in  satellite  SST  data.  In  addition,  limited  amounts  of  altimetry,  BT,  and 
other  data  are  sometimes  used.  These  data  are  sometimes  adequate  for  mapping  the  front  at  the 
surface,  but  there  are  rarely  enough  subsurface  data  for  mapping  the  three-dimensional  thermal 
structure.  The  surface  information  in  the  bogus  message  is  assimilated  into  the  thermal  analysis  at 
depth  by  means  of  feature  models  (Robinson  and  Walstad  1987;  Bennett  et  al.  1989). 

The  ring  feature  model  of  Bennett  et  al.  (1989)  is  an  example  of  a  feature  model  that  can  be 
tuned  using  ring  parameters.  This  model  initially  uses  a  typical  ring  rotation  rate,  the  observed  ring 
radius  from  the  bogus,  and  simple  dynamical  arguments  to  estimate  the  field  of  dynamic  height 
anomaly  of  the  ring.  The  thermal  structure  of  the  ring  can  then  be  modeled  by  inferring  from  the 
dynamic  height  anomaly  a  temperature  profile  at  each  point  on  a  close-mesh  grid  centered  on 
the  ring.  The  inferred  profiles  are  often  referred  to  as  synthetic  profiles.  Bennett  et  al.  (1989)  used 
an  iterative  procedure  to  infer  the  synthetic  profile.  An  alternative  approach  for  inferring  synthetic 
profiles  of  temperature,  salinity,  and  sound  speed  has  been  evaluated  by  Carnes  et  al.  (1990).  They 
derived  the  empirical  orthogonal  functions  (EOFs)  of  the  variability  of  temperature  and  salinity 
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profiles  from  a  domain  of  the  ocean,  and  then  determined  a  regression  between  the  dynamic  height 
anomaly  and  the  amplitude  of  the  EOFs.  It  is  then  straightforward  to  infer  the  EOF  amplitudes  from 
an  estimate  of  dynamic  height  anomaly  and  to  construct  a  synthetic  profile. 

Experience  has  shown  that  the  ring  rotation  rate  is  a  key  parameter  for  tuning  the  ring  feature 
model.  For  example,  consider  a  cold-core  ring  that  is  described  by  the  conceptual  model  presented 
in  Joyce  (1984)  and  Olson  (1980)  and  discussed  below.  Assume  the  radius  of  the  ring  core  is  75  km 
and  the  dynamic  height  anomaly  at  the  edge  of  the  core  is  1.5  m,  which  is  approximately  that  of 
the  ambient  water.  For  rotation  rates  of  2.0*10”^  s“’  and  2.5*10”^  s"’,  which  is  the  approximate 
range  of  measured  rates  for  cold-core  rings,  the  dynamic  height  anomaly  at  the  ring  center  is 
estimated  via  Eqs.  (1)  and  (2)  (presented  below)  to  be  0.907  m  and  0.723  m,  respectively.  The 
dynamic  height  anomaly  and,  hence,  the  subsurface  thermal  structure  are  seen  to  be  sensitive  to 
the  rotation  rate. 

Glenn  et  al.  (1990)  describe  techniques  for  estimating  several  ring  parameters,  including  the 
rotation  rate,  by  fitting  feature  models  to  satellite  SST  data,  BT  data,  and  drifting  buoy  data.  They 
found  that  a  feature  model  analysis  of  drifting  buoy  data  is  an  effective  means  for  monitoring  the 
evolution  of  Gulf  Stream  rings.  In  addition,  they  used  these  techniques  to  provide  real-time  support 
for  deep-water  drilling  operations  off  the  East  Coast  of  the  U.S. 

A  number  of  techniques  exist  for  assimilating  the  information  in  the  bogus  into  an  analysis 
via  feature  models.  In  the  implementation  known  as  OTIS  1.1  (Bennett  et  al.  1989),  the  bogus  and 
feature  models  are  used  to  construct  the  first-guess  field,  which  resembles  a  climatology  but  has 
typical  frontal  and  ring  thermal  structures  built  into  it  at  the  locations  indicated  in  the  bogus.  The 
BT  and  SST  data  are  then  assimilated  via  optimum  interpolation  into  this  first-guess  field.  In  data- 
sparse  areas,  the  analysis  will  resemble  the  first-guess  field,  which  has  schematic  descriptions  of 
mesoscale  features  at  the  appropriate  locations.  In  the  OTIS  3.0  implementation,  the  first-guess  field 
is  a  climatology.  The  feature  model  is  then  expressed  as  a  number  of  synthetic  profiles  of  tempera¬ 
ture  that  are  assimilated  along  with  BT  and  other  data  into  the  climatology  via  optimum  interpolation. 


APPROACH 

This  report  presents  a  technique  for  estimating  the  radius  and  rotation  rate  of  a  ring  from  a 
track  of  altimetry  data.  The  technique  is  applicable  even  if  the  track  does  not  cross  the  center  of 
the  ring.  In  addition,  it  is  shown  how  to  composite  estimates  from  several  such  tracks  into  a  best 
estimate  of  these  ring  parameters. 

Consider  a  circular  ring.  Observations  of  warm-core  (Joyce  1984)  and  cold-core  (Olson  1980) 
rings  indicate  that  near  the  surface  a  ring  consists  of  a  disk  in  solid  body  rotation  surrounded  by 
a  transition  zone  within  which  the  azimuthal  velocity  rapidly  decreases  to  zero.  The  azimuthal 
velocity,  v,  satisfies  the  equation 


fv  + 


V 

r 


30  _  dz 
dr  ^  dr 


(1) 


where  /  is  the  Coriolis  parameter,  r  is  the  radial  distance,  and  O  is  the  geopotential  anomaly 
between  the  sea  surface  and  a  reference  level.  The  height  (dynamic  topography)  associated  with  the 
geopotential  anomaly  is  z  =  O/g,  w  ‘re  g  is  the  acceleration  of  gravity.  Substituting  the  condition 
of  solid  body  rotation,  v  =  ht,  whe  .  w  is  the  rotation  rate  of  the  ring,  into  Eq.  (1)  and  integrating 
from  the  center  of  the  ring  outward  gives 
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Consider  a  track  of  altimeter-derived  dynamic  topography  data  that  crosses  a  disk  of  radius  R  and 
comes  within  a  distance  L  of  the  ring  center  (see  Fig.  1).  In  this  report,  c  is  used  to  represent  the 
dynamic  topography.  Denote  the  minimum  (maximum)  dynamic  topography  c  along  the  track  as  it 
crosses  the  ring  as  an  origin  and  let  /  be  the  distance  along  the  track  from  this  origin.  Beginning 
with  Eq.  (2),  the  dynamic  topography  along  the  track  can  be  expressed  as 


VV(H  +/)  2  H(M  +/)  2 

^'-^0+  2g  ^  2g  ' 


(3) 


or,  more  concisely,  as 


where 


n(vf  +f)  ^2 

2g 


(3a) 


(3b) 


and 
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Fig.  1  —  The  ground  track  of  the  satellite  altimeter  crosses  a 
ring  at  a  distance  L  from  the  center  c  of  the  ring.  The  dynamic 
topography  of  a  cold-core  (warm-core)  ring  has  a  minimum 
(maximum)  along  the  track  at  location  0,  which  is  labeled  the 
origin  in  this  context.  The  point  p  along  the  track  is  a  distance 
/  from  point  0  and  a  distance  r  from  the  ring  center  at  c.  Note 

r={r- -*■  L-)'’'. 


The  parameter  z'q  is  an  estimate  of  the 
minimum  (maximum)  dynamic  topography 
along  the  track  as  the  altimeter  crosses  the 
cold-core  (warm-core)  ring.  The  curvature 
(rate  of  change  of  slope)  of  the  dynamic 
topography  is  2*c. 

The  application  of  the  approach  begins 
by  estimating  the  dynamic  topography  along 
the  altimeter  track.  The  location  of  the  origin 
(the  dynamic  topography  minimum  or  maxi¬ 
mum  along  the  track  as  it  crosses  the  ring) 
is  then  determined  by  inspection  and  the 
dynamic  topography  is  plotted  as  a  function 
of  distance  /  from  the  origin.  Least  squares 
are  used  to  determine  the  quadratic  equation. 
Eq.  (3a),  that  best  fits  the  dynamic  topography 
as  a  function  of  /.  The  rate  of  rotation  w  is 
then  readily  estimated  from  c  and  the  known 
latitude  via 

...  - 

2 


(4) 
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The  positive  (negative)  root  is  used  in  Eq.  (4)  for  cold  (warm)  rings.  Note  that  the  estimate  of  the 
ring  rotation  rate  does  not  depend  on  the  distance  from  the  track  to  the  ring  center  as  long  as 
the  track  crosses  the  core  that  is  in  solid  body  rotation. 

The  rotation  rates  estimated  from  several  tracks  of  altimetry  data  can  be  readily  assimilated  into 
a  composite  estimate.  Consider,  for  example,  repeated  passes  of  the  altimeter  along  a  particular 
ground  track  as  a  cold-core  ring  drifts  across  the  ground  track.  The  composite  w  can  be  estimated 
by  averaging  those  estimated  values  of  tv  that  pass  a  gross  error  check.  The  uncertainty  in  the 
composite  u’  can  be  estimated  from  the  variabiluj  of  the  w  estimates  within  the  albeit  quite  limited 
distribution.  In  addition,  the  minimum  ;'o  that  passes  a  gross  error  check  is  the  best  estimate  of  c 
at  the  center  of  the  cold-core  ring. 


RESULTS 

Rotation  Rate  from  a  Single  Altimeter  Track 

In  this  section,  we  estiniate  the  rotation  rate  of  a  cold-core  ring  from  a  single  track  of  altimetry 
data.  Figure  2  is  a  plot  of  the  dynamic  topography  along  the  ascending  ground  track  A4  of  the 
GEOSAT  Exact  Repeat  Mission.  The  technique  for  estimating  dynamic  topography  is  a  modification 
of  the  approach  of  Carnes  (private  communication)  and  is  summarized  in  the  appendix.  These  data 
are  from  revolution  11311  of  the  GEOSAT  altimeter,  which  occurred  on  12  May  1987.  The  location 


Fig.  2  —  The  dynamic  topography  along  ground  track  A4  for  revolution  11311  of 
the  GEOSAT  altimeter.  A  cold-core  ring  is  located  near  37°  N-38°  N.  The  location 
of  the  minimum  dynamic  topography  is  the  origin. 
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of  a  cold-core  ring  and  the  choice  of  the  origin  are  also  indicated  on  Fig.  2.  By  comparing  this 
track  of  dynamic  topography  to  neighboring  tracks  and  by  examining  the  time  evolution  of  the  tracks, 
the  interpretation  of  the  data  in  Fig.  2  as  a  front  and  cold-core  ring  rather  than  as  a  front  and  a 
warm-core  ring  can  be  made. 

Figure  3  is  a  plot  of  all  dynamic  topography  estimates  that  lie  within  125  km  of  the  origin.  The 
"+”  signs  represent  data  from  points  that  lie  south  of  the  origin,  and  the  signs  represent  data 
from  points  to  the  north.  A  least-squares  quadratic  fit  of  all  of  the  data  that  lie  within  a  distance 
s  of  the  origin  is  made  for  several  values  of  s.  The  resulting  curves  are  shown  in  Fig.  4.  The  data 
that  lie  within  about  75  km  of  the  origin  fit  a  quadratic  curve  rather  well;  a  disk  in  solid  body 
rotation  describes  these  data  well.  Between  75  and  125  km,  there  is  the  transition  zone  within 
which  the  azimuthal  velocity  decreases  to  zero  with  increasing  distance  from  the  origin.  The  dynamic 
topography  and  fitted  curve  diverge  with  increasing  distance  within  the  transition  zone.  The 
dynamic  topography  at  a  distance  of  125  km  is  approaching  that  of  the  ambient  water. 

Table  1  presents,  for  several  values  of  s,  the  coefficient  of  determination  (Mack  1966)  and 
standard  deviation  of  the  data  about  the  fitted  quadratic.  The  coefficient  of  determination  is  the 
ratio  of  the  variance  explained  by  the  curve  to  the  total  variance.  The  coefficient  is  close  to  1.0. 
indicating  that  the  curve  fits  the  data  well  as  can  be  seen  by  inspection.  The  standard  deviation  of 
the  data  about  the  fitted  curve  is  only  on  the  order  of  a  few  centimeters.  The  data  are  from  a 
crossing  of  a  cold-core  ring  by  the  GEOS  AT  altimeter  during  revolution  11311. 
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Fig.  3  —  A  plot  of  dynamic  topography  (from  Fig.  2)  as  a  function  of  distance 
from  the  origin.  The  +  symbols  signify  points  south  of  the  origin  and  the  *  symbols 
signify  points  north. 
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Fig.  4  —  Quadratic  curves  fitted  in  a  least-squares  sense  to  the  dynamic  topography 
vs.  distance  data  of  Fig.  3.  Separate  fits  are  made  of  data  that  is  within  25,  50. 
75,  and  100  km  of  the  origin. 


Table  1  —  The  Coefficient  of  Determination  and  the  Standard 
Deviation  of  the  Dynamic  Topography  vs.  Distance  Data  About  the 
Fitted  Curve 


s  (km) 

25 

50 

75 

100 

Coefficient  of  Determination 

0.95 

0.95 

0.91 

Standard  Deviation  (cm) 

1.3 

2.8 

5.5 

As  the  data  to  be  fitted  are  restricted  to  lie  closer  to  the  origin,  the  curvature  of  the  fitted  curve 
increases.  Table  2  presents  the  parameters  Zq  and  c,  as  well  as  the  rotation  rate  w  for  each  of  the 
curves  fitted  to  the  data  and  shown  in  Fig.  4.  The  estimated  w  is  1.7*10“^  s“'  when  a  fit  is  made 
to  all  of  the  data  that  lie  within  75  km  of  the  origin.  If  the  origin  were  to  lie  at  the  center  of  the 
ring,  then  the  ring  radius  would  be  about  75  km.  The  data  are  from  a  crossing  of  a  cold-core  ring 
by  the  GEOS  AT  altimeter  during  revolution  11311. 

Repeated  Altimeter  Tracks 

In  this  section  we  consider  the  case  of  repeated  passes  of  the  satellite  altimeter  along  a  single 
ground  track  as  a  cold-core  ring  drifts  across  the  ground  track.  Figure  5  presents  the  dynamic 
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Table  2  —  Parameters  ;'o  and  c  and  the  Estimated 
Ring  Rotation  Rate  vv  for  Several  Least-Squares 
Quadratic  (;'o  +  r  /  *)  Fits  of  the  Dynamic  Topography 
vs.  Distance  Data 


topography  for  revolutions  1082.^.  1 1067.  and 
11311  along  the  GEOSAT  ascending  ground 
track  A4.  The  tracks  of  dynamic  topography 
are  offset  by  0.2  m  for  clanty.  These  tracks  are 
17  days  apart  and  are  dated  8  April  1987, 
25  April  1987.  and  12  May  1987.  The  data  for 
revolution  11311  were  considered  in  Sec  4  1. 

Figures  6  and  7  are  plots  of  the  dynamic 
topography  vs.  distance  from  the  origin  for 
tracks  10823  and  1  1067,  respecti%  els .  Also 
shown  on  each  figure  is  the  quadratic  curve 
fitted  in  a  least-squares  sense  to  all  the  data 
that  lie  within  75  km  of  the  origin.  The  curve 
is  shown  for  the  case  of  fitting  all  data  within 
75  km  of  the  origin  because  the  coefficient  of  determination  is  larger  at  s  =  75  km  than  at  .v  =  25. 
50,  or  100  km,  although  the  difference  was  minimal  for  revolution  10823.  In  the  case  of  track 
11067,  the  data  from  north  and  south  of  the  origin  did  not  agree  with  each  other  as  well  as  did 
similar  comparisons  from  the  other  two  tracks.  The  dynamic  topography  at  the  origin  is  about 
0.18  m  higher  in  the  case  of  track  10823  than  track  11067,  suggesting  that  the  track  of  altimetry 
data  crosses  the  ring  at  a  longer  distance  from  the  ring  center  in  the  case  of  track  10823.  Both  tracks 


.V  (km) 

;'()  (m) 

c  (10  '*m  km  ■) 

w  (10^  s-') 

25 

1.0 

1  247 

2.3 

50 

1.0 

i.055 

2.0 

75 

1.0 

0.870 

1.7 

100 

1.0 

0.666 

1.3 

Fig.  5  —  The  dynamic  topography  along  ground  track  A4  for  revolutions  1082.^, 
1 1067,  and  1 131 1  of  the  GEOSAT  altimeter.  The  tracks  are  offset  by  0.2  m  for 
clarity.  These  revolutions  are  consecutive  passes  along  track  A4  during  the  GEOSAT 
Exact  Repeat  Mission.  (Note  the  cold-core  ring  drifting  across  the  track  near 
36°N-38°N) 
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DISTANCE  (km) 

Fig  6  —  A  plot  of  dynamic  topography  from  GEOSAT  revolution  10823  as  a 
function  of  distance  from  the  origin.  The  +  symbols  signify  points  south  of  the 
origin  aud  the  *  symbols  signify  points  north.  Also  shown  is  a  least-squares 
quadratic  Tit  of  the  data  that  lie  within  75  km  of  the  origin. 


10823  and  11067  have  a  “bump”  in  the  dynamic  topography  relative  to  the  quadratic  curve  at 
distances  of  30-40  km. 

The  ring  rotation  rates  estimated  from  these  curves  are  tabulated  in  Table  3.  The  parameters 
are  obtained  by  making  a  least-squares  quadratic  +  of  dynamic  topography  vs.  distance 

data.  Each  fit  is  made  of  the  data  that  lie  within  75  km  of  the  origin.  Composite  1  is  a  composite 
of  all  three  tracks  while  composite  2  is  a  composite  of  tracks  1 1067  and  11311  only.  A  value  of 
j  =  75  km  is  also  used  for  revolution  11311.  The  coefficient  of  determination  for  this  case  (0.95) 
was  comparable  to,  although  smaller,  than  the  largest  value  (0.98).  However,  the  s.  bjective  judgment 
was  made  that  a  value  of  s  that  is  long  and  representative  of  the  range  of  distances  over  which  the 
conceptual  model  is  considered  valid,  is  better  than  a  very  short  value  of  s.  These  data  suggest  that 
a  single  track  of  altimetry  is  insufficient  for  making  a  confident  estimate  of  the  ring  rotation  rate. 
The  rates  estimated  from  tracks  11067  and  11311  agree  rather  well  with  each  other  but  the 
rate  determined  from  track  10823  is  only  60%  of  that  from  track  11311. 

The  most  confident  estimate  of  the  ring  parameters  is  obtained  by  compojiting  estimates  from 
several  tracks.  The  composite  w  is  the  average  of  all  of  the  values  that  pass  a  gross  error  check. 
In  this  case,  the  estimates  of  w  from  tracks  11067  and  11311  agree  well  with  each  other,  but  the 
estimate  from  track  10823  is  somewhat  of  an  outlying  estimate.  All  three  tracks  are  included  in 
composite  1  to  give  a  composite  w  of  1.4*10"^  s"*.  If  only  tracks  1 1067  and  11311  were  included 
(composite  2),  then  the  composite  rotation  rate  is  1.6*10“^  s”'.  In  addition,  all  three  values  of  cq 
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Fig.  7  —  A  plot  of  dynamic  topography  from  GEOSAT  revolution  1 1067  as  a 
function  of  distance  from  the  origin.  The  +  symbols  signify  points  south  of  the 
origin  and  the  *  symbols  signify  points  north.  Also  shown  is  a  least-squares 
quadratic  fit  of  the  data  that  lie  within  75  km  of  the  origin. 


Table  3  —  Parameters  Zq  and  c  and  the  Estimated  Ring 
Rotation  Rate  w  for  Each  of  Three  Tracks  of  GEOSAT 
Altimetry  Data  that  Crosses  a  Particular  Cold-Core  Ring 


Rev. 

z'o  (m) 

c  (10  ‘'m  km 

w  (10-5  s-l) 

10823 

1.1 

0.489 

11067 

0.766 

1.5 

11311 

0.870 

1.7 

Composite  1 

0.9 

1.4 

Composite  2 

0.9 

1.6 

are  plausible  and  the  minimum  value,  0.9  m,  is  chosen  as  the  best  estimate  of  the  dynamic  topography 
at  the  center  of  the  ring. 

As  the  ring  drifts  past  track  A4  and  across  track  A3,  estimates  of  the  rotation  rate  obtained 
using  the  data  from  along  track  A3  can  also  be  assimilated  into  the  composite  estimate  of  the 
rotation  rate.  If  the  ring  were  large  enough  to  extend  across  two  neighboring  tracks,  then  estimates 
obtained  from  both  tracks  could  be  used  in  the  composite. 
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DISCUSSION 

While  there  are  no  concurrent  altimeter  and  ring  rotation  rate  data  to  verify  these  estimates  of 
the  ring  rotatior  ate.  these  estimates  can  be  compared  to  the  measured  rotation  rates  of  other  rings. 
The  comparison  shows  that  these  estimates  are  realistic  but  perhaps  a  little  slow.  Vastano  et  al. 
(1980)  gives  the  azimuthal  speed  of  cold-core  ring  Bob  (Olson  1980)  as  2.26  km  h“*  at  a  radius 
of  30  km,  which  gives  a  rotation  rate  of  2.1*10"*'  s”’.  Richardson  ( 1983)  has  reviewed  Gulf  Stream 
rings,  including  their  rotation  rates.  He  states  that  young  rings  have  surface  swirl  velocities  of 
about  1.5  m  s“’  at  a  distance  of  30-60  km  from  the  ring  center.  The  corresponding  rotation  rate 
is  2.5*10“^  s“'  to  5.0*10"^  s"'.  Flierl  (1979)  describes  a  simple  model  of  warm-  and  cold-core 
rings.  His  model  assumes  an  inviscid,  uniformly  rotating  two-layer  fluid.  The  lower  layer  is  at  rest 
and  the  upper  layer  has  constant  potential  vorticity.  Using  the  Flierl  model  (Flierl  1979,  Fig.  10) 
gives  an  estimated  rotation  rate  of  2.4*10"^  s~*  for  a  cold-core  ring  with  a  radius  of  75  km. 

Individual  estimates  of  the  rotation  rate  can  differ  for  a  number  of  reasons.  The  geoid  could 
have  errors  in  it.  The  ring  could  be  elliptical  rather  than  circular.  The  wet  tropospheric  (water 
vapor)  correction  to  altimeter-derived  ocean  topography  can  be  large  enough  to  mimic  or  mask  true 
ocean  features  (Phoebus  and  Hawkins  1990).  Without  an  adequate  estimate  of  this  correction, 
consecutive  altimeter-derived  estimates  of  the  dynamic  topography  of  a  ring  could  show  significant 
variation  due  just  to  atmospheric  water  vapor.  These  data  did  not  have  a  wet  tropospheric  correc¬ 
tion.  The  rotation  rate  could  also  change  over  time.  The  decay  time  scale  for  available  potential 
energy,  kinetic  energy,  and  total  energy  for  an  isolated  Gulf  Stream  ring  is  on  the  order  of  500  days. 
However,  interactions  with  the  bottom  topography  and  the  Gulf  Stream  can  decrease  the  time  scale 
by  a  factor  of  10  (Olson  et  al.  1985).  Hence,  estimates  of  the  ring  rotation  rate  before  and  after  such 
interactions  can  be  quite  different. 

Estimates  of  the  ring  rotation  rate  from  other  data  sources  can  also  be  assimilated  into  the 
composite  estimate.  In  particular,  the  spiral  path  of  a  drifting  buoy  that  is  imbedded  in  a  ring  can 
be  used  to  infer  the  location  of  the  ring  center,  the  translation  speed  and  direction,  and  the  rate  of 
rotation  (see,  for  example,  Armi  et  al.  (1989)  and  Richardson  et  al.  (1989)).  A  30-day  record 
of  drifting  buoy  data  could  be  subdivided  into  three  10-day  records  and  a  rotation  rate  estimated 
from  each  10-day  period.  At  least  one  complete  ring  rotation  would  likely  occur  within  the  10-day 
period.  The  rotation  period  of  a  new  ring  is  about  2  days  and  that  for  an  aged,  less  energetic  ring 
is  about  5-10  days  (Richardson  1983).  The  three  estimates  of  ring  rotation  rate  inferred  from  the 
altimeter  could  be  averaged  with  the  three  estimates  from  the  drifting  buoy  data.  The  relative 
weights  of  the  two  types  of  estimates  could  be  determined  from  the  relative  spread  of  the  individual 
estimates  about  the  average  for  their  data  type. 

Satellite  altimetry,  satellite  IR  imagery,  and  drifting  buoys  could  complement  each  other  to 
become  key  elements  of  a  cost-effective  system  for  monitoring  ocean  rings,  especially  submerged 
culd  rings  that  have  no  surface  signature.  The  altimeter  data  could  be  used  to  determine  that  a  ring 
exists  in  an  area  and,  using  the  present  technique,  estimate  its  rotation  rate.  However,  the  distance 
of  the  altimeter  track  from  the  ring  center  is  unknown,  and  it  is  not  known  on  which  side  of  the 
altimetry  track  the  ring  center  lies.  Satellite  IR  imagery  could  be  used  in  some  cases  to  locate 
the  ring  center,  to  estimate  the  ring  radius,  and  to  determine  how  far  from  and  on  which  side  of  the 
track  the  ring  center  lies.  A  drifting  buoy  can  be  used  to  locate  the  ring  center  and  estimate 
the  rotation  rate  for  a  period  of  up  to  several  months,  including  those  periods  when  the  ring  lies 
between  the  altimetry  tracks.  However,  drifting  buoys  are  an  inefficient  means  to  initially  locate 
rings.  Used  together,  data  from  altimetry  and  drifting  buoys  can  overcome  these  shortcomings.  The 
altimeter  data  could  be  used  to  target  the  deployment  of  the  buoys.  Both  altimetry  and  drifting  buoy 
data  could  be  used  to  estimate  the  ring  rotation  rate.  Finally,  data  from  the  drifting  buoy  could  be 
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used  to  locate  the  center  of  the  ring.  The  location  of  the  ring  center  and  the  ring  rotation  rate  and 
radius  are  the  key  inputs  for  a  feature  model  description  of  the  ring. 


CONCLUSIONS 

This  report  presented  a  technique  for  estimating  the  rotation  rate  of  a  mesoscale  ocean  ring  from 
satellite  altimetry.  The  technique  was  applied  to  the  case  of  repealed  passes  of  the  satellite  altimeter 
as  a  cold-core  ring  drifts  across  the  ground  track.  The  following  are  conclusions  of  the  work: 

(1  )  A  least-squares  quadratic  curve  fits  the  altimeter-estimated  dynamic  topography  well  within 
the  core  of  the  ring. 

(2)  A  plausible  estimate  of  the  ring  rotation  rale  can  be  readily  obtained  from  the  fitted 
quadratic  curve. 

(3)  The  most  confident  estimate  of  the  rotation  rate  is  obtained  not  from  a  single  data  track  but 
from  a  composite  of  several  tracks. 


RECOMMENDATIONS 

The  following  recommendations  are  made. 

(1)  Study  the  proposed  technique  for  estimating  the  ring  rotation  rate  more  thoroughly  by 
considering  a  number  of  cases  for  which  the  rotation  rate  is  known  from  in  situ  data. 

(2)  Study  and  demonstrate  a  technique  for  assimilating  rotation  rate  estimates  from  a  satellite 
altimeter  and  a  drifting  buoy  into  a  composite  estimate. 
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APPENDIX 


This  appendix  documents  how  the  dynamic  topography  was  obtained  for  use  in  this  report.  The 
terminology  follows  that  of  Maul  et  al.  (1988). 

The  sea  surface  height  ssh  (the  distance  from  the  reference  ellipsoid  to  the  sea  surface)  is 

ssh  =  Zs-Za,  (A  1 ) 

where  is  the  height  of  the  satellite  above  the  reference  ellipsoid  as  determined  by  tracking  and 
Za  is  the  altimeter-derived  estimate  of  the  distance  from  the  altimeter  to  the  sea  surface.  It  can  be 
written  as  a  sum 


ssh  =  Zg  +  Zbc  +  Zfir  +  Ze .  (A2) 

where  Zg  is  the  geoid,  Zbc  **  ih®  dynamic  topography,  zbt  is  the  barotropic  part  of  the  ssh,  and  Zg 
is  a  sum  of  various  errors  such  as  orbit  error.  From  Eq.  (A2),  the  dynamic  topography  is  then 

ZBC  =  ^S^-h~^BT-^e-  (A3) 

The  geoid  Zg  obtained  by  Carnes  (private  communication)  is  used  as  Zg.  A  geoid  has  been 
estimated  by  Glenn  et  al.  (1991)  using  a  similar  technique.  The  Carnes  et  al.  geoid  is  computed  as 
follows:  Since  Zg  is  constant  in  time,  at  least  on  the  time  scales  considered  here,  it  may  be  estimated 
by  subtracting  a  1-year  mean  dynamic  height  anomaly,  Zgc*  the  mean  sea  surface  height  ssh 
for  the  same  period. 


Zg  =  ssh-ZBc-  (4) 

The  approximations  used  to  obtain  Fq.  (4)  are  discussed  by  Maul  et  al.  (1988).  The  first  step 
in  estimating  Zbc  is  to  use  a  water-mass- z-aoed  climatology  and  feature  models  to  construct  a  three- 
dimensional  thermal  analysis  of  the  Gulf  Stream  region  for  each  map  in  a  year-long  series  of 
weekly,  NAVOCEANO-prepared  front  and  eddy  maps.  The  dynamic  height  anomaly  was  then 
estimated  for  each  themal  analysis  and  the  ensemble  of  dynamic  height  anomaly  fields  subse¬ 
quently  averaged.  The  ssh  is  obtained  by  making  an  orbit  error  correction  to  each  track  of  sea  level 
data  and  averaging  the  ensemble  of  ssh's.  If  the  track  length  is  much  less  than  the  length  of  an 
orbit,  e.g.,  1000  km,  then  the  orbit  error  can  be  removed  by  making  a  linear  correction.  (This  geoid 
is  determined  only  along  the  tracks  of  altimetry  data  and  is  not  a  two-dimensional  field.) 

Maul  et  al.  (1988)  have  proposed  using  in  situ  inverted  echo  sounder/pressure  ^auge  and/or 
island  tide  gauge  data  to  eliminate  the  orbit  error  in  the  dynamic  topography  along  the  track.  These 
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data  would  provide  an  estimate  of  the  dynanaic  topography  along  the  track.  By  differencing  the 
uncorrected  altimeter-derived  estimate  and  the  in  situ  estimate,  the  orbit  error  can  be  explicitly 
estimated  and,  hence,  corrected. 

The  technique  used  here  to  correct  for  orbit  error,  which  was  developed  independently,  has 
similarities  with  that  of  Maul  et  al.  (1988)  and  also  Carnes  et  al.  (in  preparation)  but  differs  in 
several  respects.  In  this  technique,  a  reference  dynamic  height  anomaly  profile  is  constructed  along 
the  ground  track.  The  purpose  of  this  construction  is  to  build  a  quasi-realistic  profile  that  has  a 
sharp  gradient  at  the  location  of  the  instantaneous  front  and  resembles  climatology  elsewhere.  This 
series  is  a  reasonable  first  guess  of  actual  dynamic  height  anomaly.  The  difference  between  the 
reference  dynamic  topography  and  the  first  estimate  is  then  an  estimate  of  the  orbit  error.  A  linear 
correction  in  bias  ^  and  tilt  y  is  applied  to  z  such  that  a  least-squares  linear  fit  of  c  agrees  with  that 
of  the  reference  series.  The  distance  along  the  track  is  x.  As  a  result  of  this  correction,  the  bias  and 
tilt  of  the  corrected  z  is  the  same  as  that  of  the  quasi-realistic  estimate.  Admittedly,  this  technique 
cannot  account  for  large-scale  fluctuations  in  the  sea  surface  topography  due  to  large-scale  changes 
in  the  oceanography.  The  dynamic  topography  is  then 

Zgc=SSh-Zg-ZBT-Ze-(^  +  1x).  (5) 


In  practice,  the  baratropic  term  zbt 's  neglected. 

Figure  Al  shows  the  method  of  construction.  The  solid  line  in  Fig.  Al  is  the  1-year  mean 
dynamic  height  anomaly  obtained  by  Carnes  et  al.  (in  preparation),  which  has  a  broad  gradient 
associated  with  the  meandering  of  the  sharp  Gulf  Stream  front.  This  mean  dynamic  height  anomaly 
is  modified  by  inserting  into  the  series  a  sharp  linear  gradient  in  the  dynamic  height  anomaly  at  the 


Fig.  Al  —  The  reference  dynamic  topog¬ 
raphy  constructed  for  use  in  the  orbit 
error  correction.  The  curves  shown  are 
the  1-yr  mean  dynamic  height  anomaly 
(solid  curve),  the  mean  dynamic  height 
anomaly  with  a  linear  front  inserted  at 
the  location  of  the  instantaneous  front  (long 
dash  curve),  and  the  mean  dynamic  height 
anomaly  with  the  constructed  front  and  a 
steric  correction  (short  dash  curve). 
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observed  instantaneous  location  of  the  Gulf  Stream  front  (see  the  long  dashed  line  of  Fig.  Al).  To 
eliminate  the  climatological  gradient  in  the  dynamic  height  anomaly  in  the  vicinity  of  the  front,  the 
dynamic  height  anomaly  between  the  instantaneous  front  and  the  climatological  maximum  south  of 
the  front  is  made  a  constant  (see  Fig.  Al).  A  constant  value  of  0.85  m  is  assumed  north  of  the  front. 
Next,  a  steric  correction  to  the  dynamic  height  anomaly  is  made  (see  the  short  dashed  line  in 
Fig.  Al).  The  magnitude  of  the  steric  correction  is  determined  by  interpolating  in  time  the  magni¬ 
tude  of  the  monthly  mean  steric  change  for  Sargasso  and  Slope  waters.  The  correction  is  applied 
separately  to  the  Sargasso  and  Slope  waters.  Rings  are  not  considered  in  this  implementation  but 
could  be  readily  incorporated  into  the  algorithm. 


